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Abstract — In this contribution we present a decentral- 
ized power allocation algorithm for the uplink interleave 
division multiple access (IDMA) channel. Within the pro- 
posed optimal strategy for power allocation, each user 
aims at selfishly maximizing its own utility function. An 
iterative chip-by-chip (CBC) decoder at the receiver and 
a rational selfish behavior of all the users according to 
a classical game-theoretical framework are the underlying 
assumptions of this work. This approach leads to a channel 
inversion policy where the optimal power level is set locally 
at each terminal based on the knowledge of its own channel 
realization, the noise level at the receiver and the number 
of active users in the network. 

I. Introduction 

Interleave division multiple access (IDMA) has been 
identified as a promising multiple access technique in 
the context of cellular networks [Q and self-organizing 
networks (e.g ad hoc networks) Q. In this domain, 
distributed power algorithms play a central practical and 
theoretical role. Nonetheless, the study of distributed 
power allocation algorithms for IDMA is still an un- 
solved problem. The centralized power allocation (CPA) 
has been tackled by several authors (see 0, JH, and 
references therein). In these works, iterative multiuser 
detection/channel decoding is assumed at the receiver 
and the base station determines the power to be trans- 
mitted by each user, according to a global optimum 
criterion and typically having complete channel state 
information (CSI). In this contribution, a novel frame- 
work to tackle the DPA problem in IDMA systems 
using tools from game theory is presented. The proposed 
approach assumes iterative multiuser detection/channel 
decoding at the receiver based on the chip-by-chip (CBC) 
algorithm Q], and relies on the signal to noise ratio 
(SINR) evolution technique Q for non-heavily loaded 
systems. 

In our approach, we assume that each user allocates 
the power by maximizing its own utility function and 
by assuming a competitive selfish and rational behavior 
of the other users. The proposed utility function for a 
given user is the ratio between a power of the user's 
goodput (probability of an error-free detected frame) 
and its own transmit power. Our decentralized approach 



requires only the knowledge of the noise power at the 
receiver, the channel gain of the user of interest, and the 
number of active users in the system. It yields a channel 
inversion policy for the power allocation. This policy 
is applied to an IDMA system with a repetition code of 
length N bits. Simulations show that the allocated power 
is substantially independent of the number of users 
in typical operation conditions for a practical system. 
Therefore, the knowledge of the number of active users 
in the system becomes irrelevant as observed in the 
centralized case in 0- 

II. System Model 

Consider the uplink of an IDMA system with K 
chip-synchronous active users. Each mobile and the 
base station are equipped with a single antenna. The 
base station performs CBC iterative multiuser decoding 
based on successive interference cancellation (SIC) or 
parallel interference cancellation (PIC) as proposed in 
fill . Denote by the information bits of the length-M' 
message to be transmitted by the k th user. An identical 
code with low rate R is applied to the messages of all the 
users. The coded information bits, referred to as chips, 
are permuted by an interleaver ttj, of length M = 
chips. Each interleaver is unique in the network, i.e. 
7Tfe / tt{, Vi / k. The base band signal r (j) sampled at 
the chip-rate at the receiver is 

K 

r (j) = VPihiXi (j) + n (j) , j = . . . M - 1 (1) 

i=l 

where j is the discrete time index for the chip interval 
and hi and pi represent the channel realization and 
the transmit power of the i th user, respectively. Here, 
Xi (j) represents the chip transmitted by the i th user at 
chip interval j and n(j) is the additive white Gaussian 
noise (AWGN) process with zero mean and variance 
a 2 . For the sake of simplicity, we introduce our results 
assuming an antipodal modulation, i.e x\ (j) G {—1, 1}, 
real valued channel realizations, and real noise samples. 
The extension to more complex modulation schemes is 
straightforward. The receiver is made of an elementary 
multiuser signal estimator (ESE) based on the principle 



of maximum a posteriori (MAP) detection, a set of 
interleavers ttj. and de-interleavers 7iT , and a set of 
K single user a posteriori probability decoders (DEC). 
The decoders feed back soft information to the ESE 
module in turbo configuration to iterate and improve the 
estimations. This structure is known as CBC decoder 
and is described in HI. The SINR at the input of the 
decoder can be estimated at each iteration of the CBC 
detection by means of the SINR evolution technique J4]]. 
It is shown in [31] that in the steady state, the SINR of 
user k, denoted as 7^, converges to the solution to the 
following system of equations 
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where the super-index ss stands for steady state. The 
function /(x) G [0, 1], Vx ^ represents the amount of 
multiple access interference (MAI) that is eliminated at 
each iteration of the CBC detection 0. This function 
depends on the coding scheme and can be obtained by 
Monte Carlo simulations as described in [ 1 ]. The system 
© might have several solutions. Nonetheless, it has been 
shown in [3] that in the case of \hi\ 2 pi = \hj\ 2 pj Vi 7^ j, 
the system has a unique solution given by the fix point 
equation 



r 



(if-l)|%PW(7 ss ) + ^ 



Vfc €[!,!?] (3) 



and Yk = 7 SS . with 1 < * < K - 

Additionally, it has been shown that the SINR evo- 
lution technique predicts precisely the SINR under the 
constraint of non-overloaded systems 0, 0. However, 
we have found in our simulations that this is not the case 
for heavily overloaded systems. Under this constraint, 
the SINR evolution technique does not match the SINR 
obtained at the output of the ESE module. In order to 
exclude this case, we restrict our study to non-overloaded 
systems, i.e. we consider that ^ ^ 1 always holds. 

iii. a non-cooperative power allocation 
Game 

We define the power allocation problem as a strategic 
game denoted by the triplet {S, V, (uk) keS }, where 
S = {1, . . . , K} is the set of players or active mobiles, 
V represents the set of strategies, and Uk is the utility 
function of user k. In this work, the strategy of each user 
consists of all the possible transmittable power levels 
p k G R + , \/k G S. We follow the same line as in 0, 
and [6 ] and enforce constraints on the maximum transmit 
power later. Furthermore, the utility function u k of user 
k is defined as the ratio between the positive power Sk 



of its goodput 5(7^), i.e. the probability of an error-free 
detected packet, and its transmit poweiQ 



Uk (P,h) 
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(4) 



where p = (pi,p 2 , ■ ■ ■ ,Pk) and h = {hi, h 2 , . . . , h K ) 
are the vectors of transmit power levels and channel 
realizations, respectively. Expression (@]) generalizes the 
utility function presented in and when Sk = 1, 
V/c G S. With respect to the reference utility function 
with Sk = 1 the user k will choose sj. > 1 if it is 
interested in better performance at the expenses of a 
higher transmitted power. On the contrary, if it needs 
to save power then it will set Sk < 1. In this study, 
we assume that all the users are aware of the fact that 
they are all interested on the same QoS profile, i.e. 
s x = . . . = S K = s. 

The objective of the player k is to determine the trans- 
mit power p* k that selfishly maximizes its utility function 
Uk(p*, h) under the assumption that a similar strategy is 
adopted by all the other users. Thus, the optimum power 
is obtained as solution to the following system of 
equations: 



-fcfir- ) [sg (7fc)Pfc^7fc ~ 9 ilk) 
V/c G S. 

(5) 

Here, g(jk)' represents the first derivative of 5(7^). Due 
to the nature of the channel (non-error free channel), 
the term 9 ^„*y > V/c G S is always non-zero. Thus, the 



system in could be simplified as 

d 

sg '{lk)Pk-^—lk - 9 ilk) = V/c G S. (6) 
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Furthermore, it has been shown in |7] that 

d 7ifc / 1 + Ikf ilk 



-Ik = — 



dpk ' K Pk Vl-7fe/'(7fe) 
where 7 = (71, ... , jk) and 



q k (7), V/c G S (7) 




lkf(lk)+Y k f (lk) 



{l-llf'(lk)) (1 + 7*/ (7*)) 04 + 1), 

llf in) 



A § (1 - iir (7,)) ' 

Replacing (0 in © we obtain the system of equations 

*4 hk) Ik ( } +1 o f J? k \ ) Qkh) - 9 (7*) = 0, V/c G s 
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'Further discussions on the definition of this utility function in 
Pk = can be found in |5| and [6|. 



In ©, the variables are only 7^, \/k G S. Therefore, 
the system does not depend directly on the channel 
realizations and the transmit power of each user. This 
property was also observed in the centralized approach 

0. Moreover, in the system ([8]), all the equations are 
identical and the system is invariant to variable permu- 
tations. This implies that the solutions are also identical, 

1. e. 71 = 72 = • • • = 7k = 7- Thus, the system © 
reduces to a single equation 
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We define 2(7) = sg' (7) 79(7) — g (7) and we name 
it target function. The zeros of the target function are 
candidates to be the optimal SINR 7*. The optimal SINR 
corresponds to the SINR at which the utility function 
of each user is maximized. If the set of optimal SINR 
7* = 7*, Vi G S is known, it is possible to obtain 
an expression for the optimal transmit power level from 
expressions (24) and (25) in 0, 
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since only positive power levels are meaningful. Note 
that the expression (TTOl could be also obtained from the 
fix point equation ® as unique solution. 

Interestingly, the power allocation for a given user 
depends on the optimal SINR 7*, its channel gain, and 
the number of active users. Therefore, the knowledge of 
all the other users' channel gains in the network is not 
required. 

IV. A Decentralized Power Allocation 
Algorithm 

In this section, we propose a power allocation algo- 
rithm based on the game investigated in SectionHIll under 
the constraint of non-overloaded systems, i.e. ^ ^ 1. In 
general, the target function 2(7) is not linear and the 
solution to 2(7) = requires a numerical approach. 
The algorithm we present here is based on the secant 
method. The iterative search of the optimum SINR 7* is 
initialized by choosing two values 71 and 72 such that 

• both 71 and 72 are not lower than the minimum 
SINR required 7 rmn for reliable communications at 
rate R without considering SIC or PIC. According 

2 The constraint is widely satisfied since the function /(■y) de- 
creases rapidly with 7 Qj, ID 



to the Shannon capacity law C = ^log 2 (l + 7). 
Thus, 7 min = 2 2R - 1 and j min < 71 < 72; 
• 2(72) < 2(71)- 

Furthermore, in practical systems, the users have 
power constraints. It might happen that the optimal 
power p* k for a given optimal SINR 7* exceeds the 
maximum transmittable power p max - In this case, a user 
i, Vi G S could either transmit at the maximum power 
Pi = Pmax or not transmit pi = 0. Let us denote 
as 7 P» ax 



" the SINR achieved by the 
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user before the first iteration of the CBC algorithm 
when transmitting at the maximum power p max . Indeed, 
if 7^g ax enables reliable decoding, i.e. 7^'o ax ^ 7min> 
then the user transmits at the maximum power level 
p max . In this case, the user does not reach the optimal 
7*. However, a reliable decoding is always ensured. 

If the condition 7^g ax ^ 7min is not satisfied, the k th 
user does not transmit and is considered in outage. Note 
that, the condition over the SINR 7? g a * I s a necessary but 
not sufficient condition for a user not to be decoded. In 
fact, certain users could attain an SINR 7^g ax higher than 
7min after decoding as a result of the iterative decoding. 
However, a user cannot evaluate this possibility. In fact, 
due to the incomplete knowledge available about other 
users, it can not determine if iterative decoding is able 
to sufficiently improve the initial SINR and enable a 
reliable decoding. Then, transmitting might result in a 
waste of energy and additional interference for all the 
other users. 

Therefore, the power allocation rule is 
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otherwise 



Pk Pmax 

P* k > Pmax and 7P j x > 7„ 



(12) 

Let us denote with e > the desired accuracy to deter- 
mine 7* . The power allocation algorithm is summarized 
as follows 

1) Initialization of the secant method 

71 = Imin 

while (2(71) ^ z(7min)) then 

71 = 71 + A with A G (0, |] 
end 

72 = 71 + A with A G (0, |] 

2) Iterative step of the secant method 

while |7 i+ i - 7j| > e 

7i+i = 7i - sJPkL) * M 
end 

3) End of the secant method 

1* = 7i+l- 

4) Power allocation Determine the transmit power 
according to ([T2l . 
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Figure 1 . Average Utility Function versus the SINR in decibels for 
a standard IDMA system with K — 16 users, M = 1000 bits and 
coding rate R — j^. Solid and dashed lines correspond to the SINR 
evolution and the system simulation, respectively. 



V. Performance Assessment 

The performance of the proposed power allocation 
scheme is assessed assuming antipodal modulation and 
a repetition code with different rates R = and N = 
16, 32, 64. The channel gains and the noise are real and 
Gaussian distributed with zero mean and unit variance. 

We refer to the SINR at which the utility function 
((U) is maximized as the optimal operating point 7*. An 
approximate expression for the bit error rate when a 
repetition code with rate R = is used is 

(1 - Q (V^)) lf] ~ l 

In (fT3l) . the function Q{.) is the complementary error 
function of a Gaussian random variable. It represents the 
probability of error of one chip in the case of antipodal 
modulation. Hence, the goodput can be written as: 

g{ lk ) = {l-P e ) M . (14) 

For the coding scheme considered in this contribution, 
the function /(•) is determined by Monte Carlo simu- 
lations as described in (T). The obtained results match 
those presented in and J8] perfectly. 

In Figure [Q we plot the average utility function versus 
the SINR 7 when s = 1. The solid line shows the 
theoretical performance determined by plugging (fl4l) 
in (@]), while the dashed lines are obtained with the 
actual BER at the output of the decoder. We noticed 
in all the non-overloaded cases that the utility has a 
unique maximum. The effects of the coding rate, the 
frame length, and the number of users on the optimum 



operating point are assessed in Figures |2 13 |U and 
[5j In Figure the target function Q is plotted as a 
function of the SINR 7 for three different coding rates 
while the frame length and the number of users are kept 
constant. Lower coding rates determine lower optimal 
SINRs. In Figure |3l the target function is plotted as 
a function of the SINR 7 considering three different 
frame lengths, keeping the coding rate and the number of 
users constant. In this case, longer frames lead to higher 
optimal SINR values. In Figure HJ we plot the target 
function as a function of the SINR for several exponents 
s of the generalized utility function (0]), while keeping 
the number of users, the frame length and the coding rate 
constant. When users are interested in better performance 
rather than saving power (s > 1), the optimal SINR 
increases and vice versa. In Figure [5j the target function 
is plotted for 4, 8, and 16 active users as a function of 
the SINR when the coding rate and frame length are kept 
constant. The four lines overlap completely. Therefore, 
the variation of the number of active users K has a 
negligible effect on the optimal SINR. Then, the optimal 
operating point is practically independent of the number 
of users. This effect was also observed in the centralized 
case in (51. 

In general, a convenient solution for a non-cooperative 
game is the Nash equilibrium (NE). By definition, the NE 
is a solution such that no player is interested on changing 
its strategy since no improvement could be obtained in 
its own utility while keeping the other users' utilities 
unchanged. In this case, even though the simulations 
results show that there is a unique solution to the game 
(7*), we could not prove analytically the uniqueness 
of the NE since the functions 7(7) and 5(7) depend 
on the coding scheme. However, particular cases have 
been already studied. In the case where no interference 
cancellation is performed, i.e 7(7) = 1, the proof is 
similar to the one presented in (5j with s = 1 for the 
CDMA case. 



VI. Conclusions 

We provided a novel framework for decentralized 
power allocation in IDMA systems based on a game- 
theoretic approach. In this context, each user aims at 
selfishly maximizing its own utility function assuming a 
similar behavior is adopted by all the other users. It leads 
to a channel inversion power allocation policy where the 
optimal power level could be set at each terminal based 
on the knowledge of its own channel realization, the 
noise energy at the receiver, and the number of active 
users in the network. Interestingly, we found that under 
practical operation conditions, the knowledge of the 
number of active users becomes irrelevant as observed 
in the centralized case in (3j. 
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Figure 2. Target function 2(7) versus SINR in decibels for a standard 
IDMA system with parameters K = 16 users, M = 1000 bits. The 
solid, dot-dashed and dashed lines correspond to the coding rates 
R — j-L , R = ^ anc i ^ — _L respectively. 



Figure 4. Target function 2(7) versus SINR in decibels for a standard 
IDMA system with coding rate R = i , frame length M = 1000 
bits, K — 16 users and different values for the parameter s 
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Figure 3. Target function 2(7) versus SINR in decibels for a standard 
IDMA system with parameters K — 16 users, coding rate R = 
bits. The solid, dot-dashed and dashed lines correspond to the frame 
lengths 128, 1000 and 4096 bits respectively. 
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Figure 5. Target function 2(7) versus SINR in decibels for a standard 
IDMA system with coding rate R = i , frame length M = 1000 
bits. The solid, dashed, and dot-dashed lines correspond to K — 4, 
K = 8 and K = 16 users respectively. 
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